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The electronic structure and lattice vibrational frequencies of the newly discovered superconduc- 
tors, LiFeAs and NaFeAs, are calculated within density functional theory. We show that in the 
vicinity of the Fermi energy, the density of states is dominated by contributions from Fe 3d states. 
We also calculate the electron-phonon coupling strength and show that it is too weak to account 
for the observed values of the superconducting transition temperatures. This seems to indicate that 
the iron-based superconductors are not of the conventional type. 



I. INTRODUCTION 

A new class of layered, high-Tj, superconductors 
has been recently discovered. Kamihara et al.^i^ re- 
ported a superconducting transition temperature Tc=26 
K in fluorine-doped LaOFeAs. Shortly afterwards, it 
was found that under pressure T,, increased to 43 
K.l^l Replacement of lanthanum with other rare earth 
metals gave a series of superconducting compounds 
ReOi^ajF^jFeAs, where Re = Ce, Pr, Nd, Sm, or Gd, 
w ith transi tion temperatures close to or exceeding 50 
K! 4 | 5 | 6 | 7 | 8 | 9 | Qxygen deficient samples were also sy nthe- 
sized and found to superconduct at 55 

K.lininna Hole 

doping, through the partial substitution of La with Sr, 
or Gd with Th, w as also found to yield superconducting 
compounds .fi^lSl Using high pressure techniques, it was 
possible to increase the concentration of the F-dopanlji3 
and to synthesize superconducting compounds whe re La 
is replaced by the late rare earth elements Tb and Dy.E^E^ 
The parent compound ReOFeAs is a layered compound 
consisting of a stack of alternating ReO and FeAs layers. 
Each ReO layer consists of an 0-sheet surrounded by 
two Re sheets. Similarly, each FeAs layer consists of an 
Fe-sheet surrounded by two As sheets such that each Fe 
atom is tetrahedrally coordin ated to fou r As atoms. Neu- 
tron diffraction measurement4^ ^^ l ^° l ^^b stablish that the 
Fe magnetic moments adopt a coUinear antiferromagnetic 
(c-AFM) order whereby ferromagnetic chains are coupled 
antiferromagnetically along the direction orthogonal to 
the chains. 

Superconductivity was also discovered in a second class 
of compounds containing FeAs layers, namely AFe2As2, 
where A is an alkaline earth metal. Hole doping, by 
partial replacement of A with alkali metals, results in 
superconducting compounds with Tc reaching 38 K in 
BaFe2As2 and SrFe2As2.'^ '^^ Partial substitution 

of Fe with Co was also shown to give a superconducting 
compound with Tc—22 K.^^ Similarly to the first class, in 
the parent compounds the Fe magnetic moments in this 
second class h ave a c oUinear AFM order with a spin- 
stripes pattern .I^SMMl both classes, the Fe magnetic 
moments in the parent compounds exhibit magnetic or- 
der, at low temperature, which disappears upon doping, 
making way for the emergence of superconductivity. This 



leads to the reasonable belief that strong electronic cor- 
relations are important in these systems, and that super- 
conductivity in these co mpounds is somehow connected 
to magnetic fluctuations.Bll ^^'^^'^^ l ^^ l ^^'^^'^'^ l ^^ l ^"'^^ ! ^^ ! ^^ ! 
Indeed, the electron-phonon coupling in LaOFeAs was 
estimated to be too smalP^ to give rise to superconduc- 
tivity within the conventional BCS formulation.'^ 

Recently, a third class of iron-based superconductors 
was discovered. LiFeAs and NaFeAs were fou nd to su- 
perconduct below 18 K and 9 K, respectively.l^^EZIlSllIl 
It turns out that in these two compounds, no magnetic 
order is detected at all temperatures. In some sense, 
these two compounds are important with regards to un- 
derstanding the mechanism of superconductivity in iron- 
based superconductors, the absence of spin density wave 
(SDW) transition, on the one hand, and the relatively low 
Tc in comparison with the first two classes of iron-based 
superconductors, on the other hand, make these two com- 
pounds possible candidates for being conventional BCS 
superconductors . 

Band structure calculations, using local density ap- 
proximation (LDA) within density functi onal theory 
(DFT), were recently reported for LiFeAS.'^SIII] it was 
found that LiFeAs is semi- metallic, and that the density 
of states (DOS) near the Fermi level is dominated by the 
Fe 3d states. Thus, the electronic structure of stoichio- 
metric LiFeAs is similar to that of the parent compounds 
of the first class, with a hole cylinder at the Brillouin zone 
(BZ) center, electron cylinders at the BZ corners, and an 
electronic DOS that decreases strongly with increasing 
energy in the vicinity of the Fermi energy. 

In this work we report DFT calculations of the elec- 
tronic and lattice properties of LiFeAs and NaFeAs. In 
particular, we calculate the electron-phonon coupling 
strength and show that it is too weak to account for 
the superconducting transition temperatures observed 
in these compounds. Our calculations, together with 
previous calculations^^ of the electron-phonon coupling 
strength in LaOFeAs, seem to indicate that iron-based 
superconductors are not of the conventional type. 



2 



II. METHOD 

The electronic structure calculations are carried out 
using the all-electron full-potential linear augmented 
plane wave (FP-LAPW) method as implemented in 
WIEN2K code The exchange-correlation potential was 
calculated using the generalized gradient approximation 
(GGA) as proposed by Pedrew, Burke, and Ernzer- 
hof (PBE)JSll The radii of the muffin-tin spheres for 
the various atoms were chosen so that the neighboring 
spheres almost touch each other. We set the parameter 
HMT^max='^i whcrc Ra/t is the smallest muffin-tin ra- 
dius and ¥^max is a cutoff wave vector. The valence wave 
functions inside the muffin-tin spheres are expanded in 
terms of spherical harmonics up to Imax — 10, and in 
terms of plane waves with a wave vector cutoff ^max in 
the interstitial region. The charge density is Fourier ex- 
panded up to Gmaa:=13ao~^, where ao is the Bohr radius. 
Convergence of the self-consistent field calculations is at- 
tained with a total energy convergence tolerance of 0.01 
mRy. 

The calculation of the frequencies of the vibrational 
modes and the electron-phonon coupling parameter was 
performed using ultrasoft pseudopotentials and an ex- 
pansion of the wave function of the valence electrons 
in terms of plane waves, with an energy cutoff of 30 
Rydbergs.f^ In both the electronic and lattice calcula- 
tions, the experimental values of t he low -temperature 
lattice constants and atomic position^^^HH are used. For 
both compounds, the crystal is tetragonal with space 
group P4/nmm. In LiFeAs, the lattice constants are a 
= 3.76982 1, c = 6.30693 1, whereas in NaFeAs, a = 
3.94729A, and c = 6.99112 A. 



III. RESULTS AND DISCUSSION 

Our results for the electronic structure calculations for 
LiFeAs and NaFeAs are summarized in figures [T] and [2] 
respectively, where the electronic density of states (DOS) 
is displayed. For LiFe As, ou r calculated DOS is similar 
to that reported earher.lSOlMJin both LiFeAs and NaFeAs, 
the DOS plots show some generic features that are com- 
mon to the parent compounds of the iron-based super- 
conductors: A DOS that is dominated by the Fe 3d states 
in the vicinity of the Fermi energy, with only a small con- 
tribution from the As and alkali metal states, and that 
is strongly decreasing with energy near the Fermi energy. 
It should be noted that in a unit cell of LiFeAs, for exam- 
ple, there are two Li, two Fe, and two As atoms. Thus to 
get the total atomic DOS, the values of the atomic DOS 
shown in [l] should be multiplied by 2. The total DOS 
is the sum of the total atomic DOS and the DOS in the 
interstitial region. 

Since the DOS at the Fermi energy, N(Ei^-), is ~4 
states/eV in both LiFeAs and NaFeAs, which is not very 
small, and because of the relatively lower T^ compared 
with the other iron-based superconductors, one may won- 



der whether electron-phonon coupling may lie behind the 
mechanism for superconductivity in LiFeAs and NaFeAs. 
This notion may be given more credence by the obser- 
vation that a sodium atom is about three times more 
massive than a lithium atom, so that if the attractive 
electron-electron interaction is mediated by the alkali- 
metal atomic vibrations, then this difference in the mass 
could explain the difference in the values of Tc between 
the two compounds via the well-known isotope effect. 

To test this idea, we carried out a calculation of the 
phonon dispersion curves and the electron-phonon cou- 
pling strength in these compounds. Since the crystallo- 
graphic point group in LiFeAs and NaFeAs is T)Ah, the 
vibrational modes at F, the BZ center, are decomposed 
according to the following irreducible representations 

onon 

The acoustic modes, with vanishing frequency at F, 
the BZ center, transform according to the and 
E,j irreducible representations. Excluding the acoustic 
modes, we are left with 15 modes with nonzero frequen- 
cies; among these, the symmetric ones are Raman-active, 
while the antisymmetric modes are infrared-active. The 
calculated frequencies of the Raman- and infrared-active 
modes at the F point of the BZ are given in Table [T] The 
phonon dispersion curves in LiFeAs, plotted along high 
symmetry directions in the BZ, are shown in Fig.[3j and 
the corresponding curves in NaFeAs are given in Fig. [4] 
Our results for the phonon frequencies at the BZ center 
may be checked by Raman scattering and infrared ab- 
sorption experiments, while the phonon dispersion curves 
may be checked by neutron scattering measurements. 

We calculated the electron-phonon coupfing parame- 
ter A and found it to be 0.29 and 0.27 for LiFeAs and 
NaFeAs, respectively. For conventional superconductors, 
where the attractive electron-electron interaction is me- 
diated by phonons, the transition temperature is given by 
the Allen and Dynes'^^ modified McMillan's formulcP^ 

_ uu, -1.04(1 -f A) 

1.2 ''^^^A-M*(1 + 0.62A)^ 

where ujiog is the logarithmic average phonon fre- 
quency, expressed in degrees Kelvin, and ^* is the 
Coulomb pseudopotential parameter, usually taken to be 
~0.13. Since ojiog in LiFeAs and NaFeAs is of the order 
of 100 K, the resulting value of Tc is much less than 
1 K. We conclude that the electron-phonon coupling is 
too weak to account for superconductivity in this class of 
iron-based superconductors. 

There are some puzzling questions that beset this 
third class of iron-based superconductors. In the 
parent compounds of the first two classes of iron- 
based superconductors, magnetic order is established 
at low temperatures, where the Fe magnetic mo- 
ments adopt a coUinear antiferromagnetic (c-AFM) or- 
der; this is unequivocally confirmed by both neutron 
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FIG. 1: (Color online) Density of states (DOS) in LiFeAs. Both the total and atomic DOS are shown. The Fermi energy is the 
zero energy. Near the Fermi energy, the DOS is dominated by the Fe 3d states. 



TABLE I: The calculated frequencies, in cm ^, of the Raman- and infrared (IR)-active modes in LiFeAs and NaFeAs. The 
modes are classified by the irreducible representations (irreps) according to which they transform. 

Lj(irrep) /cm^^ 



LiFeAs 


Raman 


121(E<,) 


188(AiJ 


225(Bij,) 


240(E,) 


294(Eg) 


356(Ai5) 




IR 


228 (E„) 


276 (E„) 


277(A2„) 


338(A2„) 






NaFeAs 


Raman 


110(E,) 


176(Ai,) 


187(E,) 


199(Ai,) 


218(Bi3) 


241(E3) 




IR 


170(E„) 


183(A2„) 


233(E„) 


253(A2„) 







diffraction measurements^^E^Iinili] and DFT calcula- 
tions.32 33 34 35 36 37 38 39 40 ^^le first class, it is only 
upon electron doping through the replacement of a small 
percentage of oxygen atoms with fluorine atoms, or the 
removal of a small percentage of oxygen atoms, that 
the magnetic order is suppressed and superconductivity 
emerges. In the second class of iron-based superconduc- 
tors, magnetic order is suppressed by hole doping through 
the replacement of some alkaline earth atoms with alkali 



atoms. To better understand the situation in the third 
class of iron-based superconductors, we carried out spin- 
polarized DFT calculations on stoichiometric LiFeAs and 
NaFeAs, in addition to the calculations reported above 
for the nonmagnetic phases of these compounds. We con- 
sidered, within GGA and GGA-I-U, various spin arrange- 
ments on the Fe sites. Similarly to the case of the first two 
classes, we find that the c-AFM phase, with a spin-stripes 
pattern, is indeed the lowest energy phase. Within GGA, 




FIG. 2: (Color online) Density of states (DOS) in NaFeAs. Both the total and atomic DOS are shown. The Fermi energy is 
the zero energy. Near the Fermi energy, the DOS is dominated by the Fe 3d states. 



the energy of the c-AFM phase in LiFeAs is lower than 
the AFM phase by 0.081 eV per Fe atom (eV/Fe), lower 
than the ferromagnetic (FM) phase by 0.085 eV/Fe, and 
lower than the nonmagnetic phase by 0.123 eV/Fe. For 
NaFeAs, the energy of the c-AFM phase is lower than 
the AFM phase by 0.048 eV/Fe, lower than the FM 
phase by 0.205 eV/Fe, and lower than the nonmagnetic 
phase by 0.182 eV/Fe. The differences are even greater 
within GGA+U, where onsite Coulomb repulsion is taken 
into account. Thus, according to our DFT calculations, 
stoichiometric LiFeAs and NaFeAs should be similar to 
the parent compounds of the first two classes, and they 
should not superconduct; instead, at low temperature, 
the stoichiometric compounds should display magnetic 
order. Deviations from stoichiometry, on the other hand, 
may suppress the magnetic order, making way for super- 
conductivity, just like doping does in the first two classes 
of iron-based compounds. In the case of LiFeAs, it is in- 
deed the case that the synthesized superconducting com- 
pounds were not stoichiometric, the chemical formula be- 
ing Lii_j.FeAs.'*® The situation is less clear in the case 



of NaFeAs but we believe, on the basis of our spin- 
polarized calculations, and the absence of any detectable 
magnetic order at low temperatures, that the synthe- 
sized NaFeAs samples must also be nonstoichiometric, 
and that further studies on sample characterization are 
necessary in this regards. 



IV. CONCLUSIONS 

In conclusion, we have presented the results of elec- 
tronic structure calculations on LiFeAs and NaFeAs, 
members of a new class of superconducting compounds. 
In similarity to other iron-based superconductors, the 
density of states in the vicinity of the Fermi energy is 
found to be dominated by contributions from the Fe 
3d states. We have also calculated the Raman and in- 
frared phonon frequencies at the Brillouin zone center, 
as well the phonon dispersion curves along high sym- 
metry directions in the Brillouin zone. We have evalu- 
ated the electron-phonon coupling parameter in LiFeAs 
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FIG. 3: Phonon dispersion curves in LiFeAs, plotted along high symmetry directions of the Brillouin zone. 



and NaFeAs, and found its value to be too small to ac- 
count for the observed superconducting transition tem- 
peratures in these compounds. Our results, taken to- 
gether with previous estimates of the electron-phonon 



coupling strength in LaOFeAs, seem to suggest clearly 
that iron-based superconductors are not of the conven- 
tional type, where the attractive electron-electron inter- 
action is mediated by phonons. 
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